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Abstract 

O . 

{vq , We study the equation of state for dark energy and explicitly demonstrate that the future 

crossings of the phantom divide line i^de = — 1 are the generic feature in the existing viable f(R) 

o : ^ mode , we also explore the toe _ of to cos „ ho , „ entropy . 

illustrate that the cosmological horizon entropy oscillates with time due to the oscillatory behavior 

£"■} \ of the Hubble parameter. The important cosmological consequence is that in the future, the sign 

>-' ■ of the time derivative of the Hubble parameter changes from negative to positive in these viable 

Q_c f(R) gravity models. 
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I. INTRODUCTION 



The cosmological observations such as supernovae la 111, cosmic microwave background 

ri fi ft 

radiation |2j, large scale structure |3], and weak lensing |4| have revealed that the universe 
has been undergoing an accelerating expansion since the recent "past", which is one of the 
most challenging problems in physics today. There are two representative approaches to 
explain the late time acceleration of the universe. One is the introduction of "dark energy" 
in the framework of general relativity 5[. The other is the consideration of a modified 
gravitational theory, such as f(R) gravity p-|8j. In this study, we will concentrate on the 
later approach. 

It has been commonly adopted that a viable f{R) gravity model needs to satisfy the 
following conditions: (a) positivity of the effective gravitational coupling, (b) stability 
of cosmological perturbations 9j, (c) asymptotic behavior to the standard A-Cold-Dark- 
Matter (ACDM) model in the large curvature regime, (d) stability of the late-time de 



Sitter point [lOJ, (e) constraints from the equivalence principle, and (f) solar-system con- 
straints [ll|. Several viable models have been constructed in the literature, such as the 
popular ones: (i) Hu-Sawicki 12[, (ii) Starobinsky 13|, (iii) Tsujikawa [14J], and (iv) the ex- 



ponential gravity 15Ml8l] models with the explicit forms shown in Table [B For other viable 



TABLE I: Explicit forms of f(R) in (i) Hu-Sawicki, (ii) Starobinsky, (iii) Tsujikawa, and (iv) the 
exponential gravity models. 
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models (e.g., models in Ref. 19] ) and references, see a recent review in Ref. [8]. 

Recently, the cosmological observational data 20] also seems to indicate the crossing of 
the phantom divide u>de = —1 of the equation of state for dark energy in the near "past". 



To understand such a crossing, many attempts have been made. T 



re most noticeable one 



is to use a phantom field with a negative kinetic energy term 21] . Clearly, it surfers a 



serious problem as it is not stable at the quantum level. On the other hand, the crossing 

firr 



U,U, 



of the phantom divide can also be realized in the above viable f(R) models [12l . Il7l . Il8l . l22j ] 
without violating any stability conditions. This is probably the most peculiar character o : 
the modified gravitational models. Other f(R) gravity models with realizing a crossing [q. 123 1 
as well as multiple crossings 24[ of the phantom boundary have also been examined. 

However, most of the studies in f(R) gravity have been focused on the past. In this 
paper, we would like to explore the future behaviors of the universe. In particular, we show 
that the viable f{R) models generally exhibit the crossings of the phantom divide in the 
"future" too. In addition, we investigate the future evolution of the cosmological horizon 
entropy and demonstrate that the cosmological horizon entropy oscillates with time due to 
the oscillation of the Hubble parameter. 

We use units of Jcb — c — h — 1 and denote the gravitational constant 8nG by n 2 = 
8tt/M p1 2 with the Planck mass of M P1 = G~ 1/2 = 1.2 x 10 19 GeV. We assume the flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time with the metric, 

ds 2 = -dt 2 + a 2 (t)dx 2 , (1) 

where a{t) is the scale factor. 

The paper is organized as follows. In Sec. II, we explain f(R) gravity and derive the 
gravitational field equations. We investigate the cosmological evolution in Sec. III. Finally, 
conclusions are given in Sec. IV. 

II. f{R) GRAVITY 

The action of f(R) gravity with matter is given by 

,4 n-J(R) 



I — Id Xy/—g 2 h 1 matter \9^vi T ma tter) i (^) 

where g is the determinant of the metric tensor g^, J ma tter is the action of matter which is 
assumed to be minimally coupled to gravity, i.e., the action / is written in the Jordan frame, 
and T matter denotes matter fields. Here, we use the standard metric formalism. By taking 
the variation of the action in Eq. ([2]) with respect to g^ u , one obtains [7( 

FG, U = K 2 T^ ttcr) ~ \d,u (FR -f) + V M V,F - g^DF , (3) 



where G^ u = R^ u — (1/2) g^ u R is the Einstein tensor, F(R) = df(R)/dR, V M is the covariant 
derivative operator associated with g^, D = g^V^Vj, is the covariant d'Alembertian for a 



scalar field, and T, 



(matter) 



[XV 



is the contribution to the energy-momentum tensor from all perfect 



fluids of matter. From Eq. Q, we obtain the following gravitational field equations: 



3FH 2 = K 2 p M + - {FR -f)- 3HF , 



-2FH = k 2 (p M + P M ) + F-HF 



(4) 
(5) 



where H = a/a is the Hubble parameter, the dot denotes the time derivative of d/dt, and 
Pm and Pm are the energy density and pressure of all perfect fluids of matter, respectively. 

III. COSMOLOGICAL EVOLUTION 



We analyze the cosmological evolution of dark energy including those of non-relativistic 
matter (cold dark matter and baryon) and radiation by solving Eq. (jlj) and R = 
6 ( 2H 2 + H). In particular, we examine the equation of state for dark energy, given by 

wde = Pde/pde, 





1 


Pde — 


^ 


-Pde = 


1 



- (FR - /) - 3HF + 3(1-F)H 2 



- (FR - f) + F + 2HF - (1 - F) (2H + 3H 2 ) 



(6) 



in the future, i.e., the region with the redshift z < 0. Furthermore, we investigate the 
evolution of the cosmological horizon entropy. It has been shown that it is possible to obtain 
a picture of equilibrium thermodynamics on the apparent horizon in the FLRW background 
for f(R) gravity due to a suitable redefinition of an energy momentum tensor of the "dark" 
component that respects a local energy conservation in Ref. 25[. In this picture, the horizon 
entropy S = A /(AG) 26(] where A is the area of the apparent horizon, is simply expressed 
as [25] 

We note that in the context of modified gravity theories including f(R) gravity, a horizon 
entropy S associated with a Noether charge has been proposed by Wald 27|, given by 28 ] 
S = F(R)A/ (AG). 
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FIG. 1: Future evolutions of 1 + wde as functions of the redshift z in (i) Hu-Sawicki model for 
p = 1, c\ = 2 and c% = 1, (ii) Starobinsky model for n = 2 and A = 1.5, (iii) Tsujikawa model for 
fj, = 1 and (iv) the exponential gravity model for /3 = 1.8, respectively. The thin solid lines show 
1 + u/de = (cosmological constant). 



In Figs. 1, 2 and 3, we depict the future evolutions of I+wde; H = H—Hf with H = H/H 
and Hi = if(^ = —1)/Hq, and S = S — St with 5 = S'/S'o and S'f = S'^ = — l)/So 
as functions of the redshift z = l/o — 1 in (i) Hu-Sawicki model for p = 1, ci = 2 and 
C2 = 1, (ii) Starobinsky model for n = 2 and A = 1.5, (iii) Tsujikawa model for \x = 1 
and (iv) the exponential gravity model for (3 = 1.8, respectively. Here, So = n/ (GHq) is 
the present value of the horizon entropy S, Hq is the current Hubble parameter, and the 
subscript T denotes the value at the final stage z — — 1. Note that the present time is 
z = and the future is — 1 < z < 0. The parameters used for each model in Figs. 1-3 are 
the viable ones 29|, |30|. Several remarks are as follows: (a) the qualitative results do not 
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FIG. 2: Future evolutions of H = H - H { with # = #/#(> and F f = H(z = -1)/H as functions 
of the redshift z. Legend is the same as Fig. 1. 

strongly depend on the values of the parameters in each model; (b) the evolutions of the 
Wald entropy 5* are similar to S 18|, |25[ in the models of (i)-(iv); and (c) we have studied 
the Appleby-Battye model 3l|, which is also a viable f{R) model, and we have found that 
the numerical results are similar to those in the Starobinsky model of (ii) as expected. 

We note that the present values of wde(z = 0) are -0.92, -0.97, -0.92 and -0.93 for the mod- 
els of (i)-(iv), respectively. These values satisfy the present observational constraints [33]. 
Moreover, a dimensionless quantity H 2 / ( K 2 pin /3] can be determined through the numer- 
ical calculations, where pm is the energy density of non-relativistic matter at the present 
time. If we use the observational data on the current density parameter of non-relativistic 
matter fim = Pm /p c J lt = 0.26 with p c J lt = 3H^/k 2 33|, we find that the present value of the 
Hubble parameter Hq = H(z = 0) is Hq = 71km/s/Mpc 33[ for all the models of (i)-(iv). 
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FIG. 3: Future evolutions of S = S — S{ with 5 = S/Sq and Sf = S(z = —l)/So as functions of 
the redshift z. Legend is the same as Fig. 1. 



Furthermore, (H f ,S f ) = (0.80, 1.6), (0.85, 1.4), (0.78, 1.7) and (0.81, 1.5), for the models of 
(i)-(iv), respectively, where H{ = H(z = —1)/H and S{ = S(z = — l)/5o- 

It is clear from Figs. 1-3 that in the future (—1 < z < —0.74), the crossings of the 
phantom divide are the generic feature for all the existing viable f{R) models. By writing 
the first future crossing of the phantom divide and the first sign change of H from negative 
to positive as z = z cross and z = z p , respectively, we find that (z, 



cross; ^p)a 



-0.83,-0.98) 



(-0.76,-0.82),, 
-0.79, — 0.80)jjj and (—0.74, — 0.80)j t ,, where the subscript a represents 



the ath viable model. The values of the ratio 5 = f^m/RoE at z 



and z = z p are 



Zp)) a = (5.2 x 10~ 3 ,2.1 x 10~ 3 )i, (1.7 x 10~ 3 ,4.8 x lO" 6 )^, (4.1 x 



10- 3 ,3.1 x 10"% and (6.2 x 10~ 3 ,2.8 x 10~ 3 ) M „ where Q DE = Pde/p^L and Q m = pj p^{ 



(o) 



(0) 



crit 



are the density parameters of dark energy and non-relativistic matter (cold dark matter 



and baryon), respectively. As z decreases (—1 < z < —0.90), dark energy becomes much 
more dominant over non-relativistic matter (5 = f2 m /f2DE ^ 10 -5 ). As a result, one has 
wde ~ w e ff = — 1 — 2H j (3H 2 ) = Ptot/Ptot) where w e g is the effective equation of state for 
the universe, and p tot = pde + Pm + Pr an d -Ptot = -Pde + Pt ar e the total energy density 
and pressure of the universe, respectively. Here, p m ( r ) and P r are the energy density of non- 
relativistic matter (radiation) and the pressure of radiation, respectively. The physical reason 
why the crossing of the phantom divide appears in the farther future (—1 < z < —0.90) is 
that the sign of H changes from negative to positive due to the dominance of dark energy 
over non-relativistic matter. As u>de ~ w e $ in the farther future, u>de oscillates around the 
phantom divide line wde = — 1 because the sign of H changes and consequently multiple 
crossings can be realized. We remark that since S oc H~ 2 , the oscillating behavior of S 
comes from that of H. However, it should be emphasized that although S decreases in some 
regions, the second law of thermodynamics in f(R) gravity can be always satisfied because S 
is the cosmological horizon entropy and it is not the total entropy of the universe including 
the entropy of generic matter. It has been shown that the second law of thermodynamics 
can be verified in both phantom and non-phantom phases for the same temperature of the 



universe outside and inside the apparent horizon in Ref. 32]. 

Finally, we mention that in our numerical calculations, we have taken the initial conditions 
of z = 8.0, 8.0, 3.0 and 3.5 for the models of (i)-(iv) at z — z , respectively, so that 
RF'(z = zq) ~ 10~ 13 with F' = dF/dR, to ensure that they can be all close enough to 
the ACDM model with RF' = 0. We note that in order to save the calculation time, the 
different values of zq mainly reflect the forms of the models, i.e., the power- law types of (i) 
and (ii) and the exponential ones of (iii) and (iv). It is clear that the smaller RF'(z = z ) is, 
the closer the model to ACDM. However, the cut off of RF'(z = zq) ~ 10 -13 is assumed to 
evade the divergence of the calculations in the computing program. We have also checked 
the results under the initial conditions with R, f, F, and F' as the ACDM values at the 
same redshift up to zq = 4.0 in the models of (i)-(iv) and found that they are qualitatively 
similar. We remark that at z = Zq, wde = — 1- 

Since R/R c ^> 1 in the high z regime (z ~ zq), the value of the combination ^R c is set 
as 7-R c ~ ISH^Qia , where (7, R c ) corresponds to (ci,i?Hs) 5 (A, i?s) 5 (p> -^t) an d ((3, Re) for 
(i) Hu-Sawicki (c-2 = 1), (ii) Starobinsky, (iii) Tsujikawa and (iv) the exponential gravity 
models, respectively, and f2m = Pm /p CT ( t = 0.26 33[. The reason is as follows. In the high 



z regime (z ~ zo), R/R c ^> 1, in which f(R) gravity has to be very close to the ACDM 
model, 1 R C ~ 2A = 2k 2 (pde//^) (pm/pSt) pS = 6 (pde/p^) H$n£\ where A is the 
effective cosmological constant in the limit of R/R c ^> 1. As an initial condition, we take 
(pDE/pm } ) = 3.0 llSj. Thus, we obtain 7 i? c ~ 18H$Q$ . 

Since the initial value at z = z of the variable yu = /Ode/Phi introduced in Ref. [18J 
is an arbitrary one, we have chosen it by using the observational data 33[ at the present 
time as yu (z = Zq) ~ O^/fiL — 3.0. The physical reason is as follows. By examining 
the cosmological evolutions of yn and wde as functions of the redshift z for the models, we 
have found that yu (z = 0) is close to its initial value of yu {z = Zq). This is because in the 
higher z regime, the universe is in the phantom phase (u>de < — 1) and therefore, pde and yu 
increase (since yu oc Pde) 5 whereas in the lower z regime, the universe is in the non-phantom 
(quintessence) phase (u>de > — 1) and hence they decrease. Consequently, the above two 



effects cancel out. To check our numerical results, by using the way in, e.g., Ref. [30], 
we have made another numerical calculation for the model of (iv) as follows. First, by 
taking a testing value of 7_R C at an initial value of zq, we analyze the cosmological evolution 
numerically. Next, we find the adequate value of zq so that the current value of yn {z = 0) 
can be equal to the present observational value of p^/pm = ( 1 — ^m ) /^m = 2.85, 
where p^ is the present energy density of dark energy. We also obtain the initial value of 
yjj (z = z ) and the value of R c ~ 6'~f~ 1 yH (z = Zo)m 2 with m 2 = H^Qin . As a result, we 
have found that yn [z = z ) = 2.72 with z = 3.0 for yu (z — 0) = 2.85. Furthermore, we 
have confirmed that the obtained results with the above method are qualitatively similar 
to the ones shown in Figs. 1-3. Clearly, our results are not sensitive to the initial values 
of z and yu [z = z ). We note that the cosmological evolution of %e as a function of the 



redshift z for the model of (iv) is given in Fig. 5 of Ref. 18J. On the other hand, the initial 
condition of dyn/dlna (z = z ) = is due to that the f(R) gravity models at z = Zq should 
be very close to the ACDM model, in which dyn/dlna = 0. 

IV. CONCLUSIONS 

In the present paper, we have explored the future evolution of wde in the viable f(R) 
gravity models and explicitly shown that in the future the crossings of the phantom divide 
are the generic feature in these models. We have also investigated the future evolution of 



the cosmological horizon entropy and demonstrated that the cosmological horizon entropy 
oscillates with time because the Hubble parameter also does. The new cosmological ingre- 
dient obtained in this study is that in the future the sign of H changes from negative to 
positive due to the dominance of dark energy over non- relativist ic matter. This is a common 
physical phenomena to the existing viable f(R) models and thus it is one of the peculiar 
properties of f(R) gravity models characterizing the deviation from the ACDM model. 
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